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The excited states of ^"^B were studied with a method of antisymmetrized molecular dynam- 
ics(AMD). The theoretical results suggest that the intruder states with large deformations construct 
' the rotational bands, K^ = 3/2~ and K^ = 1/2^, starting from 5 MeV and 8 MeV, respectively. 

\ The neutron structure of the K^ — 3/2~ is analogous to the intruder ground state of ^^Be. In the 

. predicted K^ = 1/2^ band, we found very exotic structure with a proton intruder configuration. 

fSJ ' This proton intruder state has a larger deformation than superdeformation. The band-head 1/2^ 

state is assigned to the 1/2"*" (4.83 MeV), which was experimentally suggested to be the proton in- 
truder state because of the strong production via the proton-transfer to the ^■^Be(O^) state in the 
*He(^^Be,"B7)X experiments. 
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I. INTRODUCTION 



In the recent progress of experimental and theoretical researches of unstable nuclei, various exotic phenomena have 
^ ' been discovered. As well known, one of the attractive subjects in neutron-rich nuclei is the breaking of neutron magic 
^ ' number such as iV = 8 and = 20 suggested in neutron-rich p-shell and sc?-shell nuclei. The breaking of the A^ = 8 
' shell in ^^Be has been well known for a long time because of the parity inversion of the ground state, J"^ — 1/2+. 
^ For ^^Be, various experimental researches have been recently achieved [l|, 0, to study the properties of the ground 
^ , band, and the breaking of A^ = 8 shell closure has been established. Concerning the mechanism of the N = 8 shell 
' breaking, the recent observation [Sj] of the significant c?-wave component in the ^^Be ground state is the direct probe for 
, the deformation, which should be one of the essential factors for the breaking in Be isotopes suggested in theoretical 
' calculations [1, H, H, 0, 11] ■ In the theoretical side, many kinds of microscopic calculations have been performed to 
[ investigate neutron-rich nuclei with the breaking shell closure. In case of neutron-rich Be isotopes, various properties 
\^ ' have been successfully described by many groups from a point of view of cluster (see references in Ref.|!8]). As a 
result, it is considered that molecular orbital structure with large deformation is a key for the breaking of A^ = 8 
magic number in neutron-rich Be. 

In the molecular orbital picture. Be isotopes are described by 2 a clusters and valence neutrons which occupy the 
molecular orbitals formed around the 2a cored, [13, [IH- The molecular orbitals are given by a hnear combination 
of p-orbitals around the a clusters, and they are associated with the orbitals in two-center shell model[l2|- In the 
^ . , spherical limit, the negative-parity 7r-orbitals are the lowest for the valence neutrons. On the other hand, with the 
■ development of the 2a clustering, so-called cr orbital, which is the longitudinal positive-parity orbital, comes down 
because its kinetic energy decreases. Finally, the inversion of the tt and a orbitals occurs in a system with well- 
developed 2a clustering. The inversion of the molecular orbitals corresponds to the parity inversion of the 'p' and 
'sd' orbitals because the tt and a orbitals in the molecular orbital model are associated with the p and sd orbitals in 
the deformed shell model, respectively. The intruder state of ^^Be is written by the configuration with two neutrons 
in the cr-orbital. If the 2a structure develops enough in ^^Be, the intruder state may become the ground state. This 
necessarily causes the significant mixture of d-wave component, which is consistent with the recent measurements 
0. Moreover, most of the low-lying states in neutron-rich Be isotopes can be well described by the molecular orbital 
structured, M. IM fiol. Hi. IS ^3, lH, [H, [I3|i and a variety of cluster states has been predicted in excited states of 
"Be and ^^Be [i,0,lllnni8[lli7 These facts indicate that cluster aspect is one of essential features in light unstable 
nuclei as well as in light stable nuclei. In particular, cluster structure is favored in neutron-rich Be, where a variety 
of exotic structure arises due to the the formation of 2 a clusters and the molecular orbitals. 
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Let us turn to structure of B isotopes, which have a proton number Z = 5 larger by one than Be isotopes. In 
contrast to the situation of ^^Be, ^^B is considered to have the normal ground state with a neutron p-shell closed 
configuration. Comparing with the intruder configuration of the ground state in ^^Be, this shows that the N = 8 magic 
number is restored in ^'^B due to the additional proton. Thus, the additional proton gives drastic structure change of 
the ground state, however, it is still natural to expect that intruder states may appear in the excited states of ^"^B. 
If a cluster structure can develop also in B as well as Be, the intruder states may lie in low excitation energy region, 
because such states can be stabilized in a similar way to the ground state of ^^Be. The intruder states may have large 
deformation and construct rotational bands. Cluster features in neutron-rich B isotopes have been systematically 
studied in a chain of B isotopes with a molecular orbital model(MO) and a method of antisymmetrized molecular 
dvnamics f AMD) [20l. [2]| . These studies have been concentrated on the structure of ground states, and have shown that 
^^B is most spherical among B isotopes, while developed cluster structure has been predicted in further neutron-rich 
B like ^^B and ^^B. These results suggest a trend of two-center cluster structure in Z = 5 systems. 

In this paper, we investigated deformed states of ^^B by performing microscopic calculations of the ground and 
excited states of ^■^B. We focused on their cluster aspect. In particular, the molecular orbital structure is our major 
interest. In the present study we applied an AMD method, which is known to be a powerful approach to investigate 
cluster structure of unstable nuclei [20l. [2ll [2^ . The present method is the same as that applied to the studies of 
^°Be,^^Be, and ^^Be[l, 0, [13 • Namely, we performed variation after spin-parity projection within the framework of 
AMDf23^. This method has been proved to successfully describe various properties of the ground and excited states 
of Be isotopes. 

The paper is organized as follows. In the next section, we briefly explain the theoretical method of the present 
work. Results and discussions are given in IIIII Finally, we give a summary in IIVI 



II. FORMULATION 



We performed e nerg y variation after spin parity projection (VAP) within the AMD model space, as was done in the 
previous studies 0, Ea HI]. The detailed formulation of the AMD method for nuclear structure study is described in 
2C ^ 2l|, [l^, ■ In particular, the formulation of the present calculations is basically the same as that described in 



An AMD wave function is a Slater determinant of Gaussian wave packets; 



^AMdCZ) = ^=A{(pi,ip2,---,fA}, (1) 

where the ith single-particle wave function is written by a product of spatial((/)), intrinsic spin(x) and isospin(T) wave 
functions as, 

= (I^XiXiTi, (2) 
<^x.(r,) cc exp{-j.(r,-^)2}, (3) 

X^ = (^+eOXT + (^ -6)Xi- (4) 

and Xi are spatial and spin functions, and is iso-spin function which is fixed to be up(proton) or down(neutron). 
We used a width parameter i/ = 0.18 fm~^, which is chosen to be the optimum value for ^■^B. Accordingly, an AMD 
wave function is expressed by a set of variational parameters, Z = {Xi, X2, • • • , Xyi, ^1, ^2, ■ • ■ , Ca}- 

For the lowest state, we varied the parameters Xj and ~ I A) to minimize the energy expectation value 
of the Hamiltonian, ($|i?|$)/($|<I>), for the spin-parity projected AMD wave function; $ = P//^,$amd(Z). Here, 
P'mk' the spin-parity projection operator. Then we obtained the optimum solution of the parameter set; Z'/'^ for 
the lowest J'^ state. The solution Z;^^ for the nth state are calculated by varying Z so as to minimize the energy 
of the orthogonal component to the lower states. 

After the VAP calculations for the states with respect to various J, n and tt = ±, we obtained the optimum 
intrinsic wave functions, $amd(Z;^'^), which approximately describe the corresponding states. After the VAP, we 
superposed the spin-parity eigen wave functions projected from all the obtained AMD wave functions. Namely, we 
determined the final wave functions for the states as, 

\Jn) = 5]c^"(i^, J.,^^,/cO|P,ti^<i>AMD(Zfr')). (5) 
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III. RESULTS 

We adopted the same effective nuclear interaction as tliat used in the study of Be isotopes [1, 0] , which consists 
of the central force, the spin-orbit force and the Coulomb force. The Majorana parameters in the MVl force are 
m = 0.65, and the strengths of the spin-orbit force are uj = —ujj = 3700 MeV. The VAP calculations of AMD 
using these interactions reproduce well the breaking of neutron magic number N ~ 8 in ^^Be and ^^Be. With this 
interaction the calculated binding energies of ^^Be and "'^^B are 61.9 MeV and 76.4 MeV, which underestimate the 
experimental values, 68.6 MeV and 84.5 MeV, respectively, however, we adopted this parametrization because the 
energy levels of the excited states in ^°Be, ^^Be and ^^Be are reasonably reproduced. 

The basis AMD wave functions were obtained by the VAP for the ground and excited states of ^^B. The number 
of the basis AMD wave functions in the present calculations are 23. The initial wave function in the energy variation 
was randomly chosen for J < 5/2 states. For J > 7/2 states, we started the variational calculation from the initial 

wave function projected from the obtained wave function $amd(Z^' ') for the J' < 5/2 states. These independent 
AMD wave functions were superposed to calculate the final wave functions. 



A. Energies and deformation 

^•^B is a nucleus with neutron magic number N = 8, and its ground state is the 3/2~ with normal configuration 
of p-shell closure. Above the 3/2~ ground state of ^^B, it is experimentally known that many states exist in the 
excitation energy Ex > 3.5 MeV region with high level density. Unfortunately, spins and parities of most of these 
states are unknown. Recently, the state at 4.83 MeV has been assigned to be a 1/2+ state by ''He(^^Be,"'^^B7)X 
experiments [26i]. Because of its strong production via proton-transfer to the ^^Be(0"'") state, this excited state is 
suggested to be a proton intruder state. 

The calculated energy levels of the negative- and positive-parity states of ^■^B are shown in Fig. [1] In addition to 
the ground 3/2|j~ state, we obtained many excited states with various J'^ in the region E^. > 4: MeV. These states may 
correspond to the observed levels in this energy region. In the excited states, we found three largely deformed bands, 
K"^ = 3/2~, = 1/2+ and K"^ = l/2~ (solid lines). These bands are composed of intruder states or well-developed 
cluster states. In particular, the band-head state 1/2^ of the = 1/2+ is the proton intruder state with a large 
deformation, and hence this should be assigned to the experimental l/2+(4.83 MeV) state. The K"^ = l/2~ band 
was obtained by the spin-parity projection and the diagonalization of the obtained wave functions, though the VAP 
calculations were not performed for the corresponding states. This band is dominantly the a-^Li cluster state. The 
intrinsic structures of these deformed states are discussed later in detail. As for other excited states (disconnected 
filled circles), intrinsic deformation of the major AMD wave function $amd(Z;{'^), which dominates the final wave 
functions | J^), is small or as large as normal deformation at most. These excited states are regarded to be dominated 
by Ohio or and neutron Ihuj configurations. 

Figure [2] shows density distribution and deformation parameters of the major AMD wave functions, ^amdC^i ), 

*amd(Z?/'"), $amd(Z2/'"), $amd(Z^/'"), $amd(z1/'+) and $amd(Z^/'+), which were obtained by the VAP for 
the corresponding states. The ground state (3/2j^) has the most spherical shape(Fig. [^fb)), due to the neutron 
p-shell closure. This is consistent with the previous work by AMD [20,]. In the 1/2^^ state(Fig. [^^a)), a three-center 
cluster core structure appears. The core clusters are an a with two valence neutrons, a triton and an a. This state is 
approximately regarded as the S'C/(3)-limit cluster state though the spatial cluster development is somehow contained. 
The similar three-center cluster structure is found also in the 5/2^ state. In the 3/2^ and the 5/2 J" states (Fig.[2Jc) 
and (d)), we found remarkably deformed structure with developed cluster cores. These states are the members of 
the = 3/2~ band which starts from Ex = 5 MeV. It is interesting that such a largely deformed band appears 
only 5 MeV above the ground state, even though this nucleus has neutron magic number iV = 8. The J"^ — 5/2^ 
state of this band is on the yrast line. Moreover, it is striking that a further large deformation arises in the 1/2'^ 
state(Fig. [5{e)), which is the band-head state of the = 1/2+ band. The deformation /3 — 0.74 of this state exceeds 
the value for superdeformation and is close to the value 0.9 for hyperdeformation. In the 7/2^ state, we obtained the 
well-developed cluster structure like ^Li-|-Q;(Fig. [2ff)). In the final wave functions after the superposition, these two 
components of the largely deformed state (Fig.[2je)) and the ^Li+a cluster state (Fig. [2jf)) constitute the rotational 
band, K"" ^ 1/2+, with a mixing of them. The \l/2l), \3/2^), |5/2j) states are dominated by P^+^$amd(Z}/^+) 
in about 90%, 65% and 60%, respectively. On the other hand, the I7/23) and |ll/2+) contains major percentage 

-f^if m'^amd(Z3^^~'')) with remarkable ^Li-t-a cluster structure, while other states in the = 1/2+ band are the 
mixture of these two components. It indicates that the weak-coupling cluster feature is enhanced in high spin region 
of the K'^ = 1/2+ band. The negative-parity states Pj[^^<I>amd(Z3^^^) projected from the ^Li+a cluster structure 
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construct the K'^ = 1/2 band. 
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FIG. 1: Excitation energies of the negative- and positive-parity states of ^^B calculated by superposition of the basis wave 
functions. Filled circles are the energies of the JZ states, for which the VAP calculations were done. Open circles are the 
energies of the JZ states, which were obtained by diagonalization of Hamiltonian by superposing wave functions, but the VAP 
calculations were not performed for the corresponding states. 



B. Cluster feature in the = 3/2" and K" = 1/2+ 

In this subsection, we discuss cluster features of the deformed bands, = 3/2~ and K"^ — 1/2+ , by analyzing 
the single-particle orbitals in the intrinsic states. Hereafter, we mainly analze the major AMD wave functions, 
*AMD(zf $AMD(Zy'+) and $amd(Z^/'+) obtained by VAP for the JZ = b/2^{K- = 3/2-), J- = l/2+(i^- = 
1/2+) and — 7/2'^{K'^ = l/2~) states, and compare them with that for the intruder ground state — 0^ of 
i2Be. 

First, we give single-particle energies in Fig. [31 In these deformed states, the level structure of the single-particle 
wave functions shows a feature of the 2a + p + An structure rather than that of spherical shell structure. That is to 
say, the lowest four proton orbitals and the lowest four neutron orbitals form the 2a core, while the higher orbitals 
correspond to one valence proton and four valence neutrons. 

Next, we illustrate the density distribution of the single-particle wave functions for the valence nucleons in Fig. 3] 
The total densities of protons and neutrons are also given as well as total matter densities in the figure. Generally 
speaking, the matter densities show two-center structures in all these deformed states. However, the behavior of 
the valence nucleons are different among these three states, = 5/2^{K'^ — 3/2"), — 1/2^{K^ — 1/2+) and 
J: = 7/2+(X- = 1/2-). 
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FIG. 2: Density distribution of the ground and excited states of ^^B. The intrinsic density of the dominant AMD wave function 
in the | J^) is shown. The density is integrated along the axis perpendicular to the plane. The deformation parameter /3[2J] of 
matter density for the intrinsic state are also given. The box size is 10 fm. 



In the J'^ — 5/2i{K'^ — 3/2^) state (Fig. [4lja)), two valence neutrons occupy an approximately positive-parity 
orbital, and the other two neutrons and a proton are in orbitals with dominant negative-parity component. Since 
the negative- and positive-parity orbitals of the valence nucleons are associated with the p-orbitals and sc?-orbitals, 
respectively, we can roughly describe the states in the K"^ = 3/2~ band by the neutron 2Uuj excited configurations. 
Let us turn to the molecular orbital features. The positive-parity orbital of the last two neutrons is largely deformed 
and has nodes along the 2a directionflongitudinal axis). This orbital well corresponds to the so-called a orbital in 
the molecular orbital picture [1, d, H, [13, [U • It has been already known that the cr-like orbital of valence neutrons 
appear in various Be isotopes (see references in [^). In the AMD study p?!, it has been revealed that the ground state 
of ^^Be is dominated by the intruder state with two neutrons in the a orbital, which reduce kinetic energy due to the 
developed 2Q;-core structure. In Fig.|3]^4), density distribution and single-particle orbitals of ^^Be(Of ) are shown. As 
seen in the figure, the last two neutrons in the ^^Be(O^) state occupy the a orbital. The point is that the neutron 
structure of the ^^B{K'^ = 3/2~) band is very similar to that of the ^^Be{Oi). Therefore, we conclude that the 
— 3/2" is the band of the intruder neutron 2hui states, and interpreted as ^^Be{Of )+p, where the ^^Be{Of) has 
the intruder configuration and the additional proton strongly couples to the deformed core. It is also interesting that 
the additional proton in the normal p-shell affects a change of the deformation of total density, which results in the 
smaller deformation of the '^^B{K'' = 3/2~) than the ^^BeiK"" = Oj'"). 

In the band-head 1/2^ state of the K"^ = 1/2+ band (Fig. UKb)), Ihui excitation occurs in the proton shell. Namely, 
the last proton occupy a tr-like orbital, which is quite similar to that of the highest neutron orbital in the = 3/2— 
band and also that in the ^^Be(O^). This is extraordinary configuration because excitations are naively expected in 
the neutron side in case of neutron-rich nuclei. It can be understoo d by the lowering mechanism of the a orbital due 
to the developed two-center structure. In the two-center shell modelflj], the energy of the a orbital comes down with 
the increase of two-center distance d. On the other hand, the negative-parity levels(7r) originating in the p orbitals 
split and four of the negative-parity levels go up as the two-center distance increases. As a result, the inversion of the 
positive-parity a orbital and the negative-parity tt orbitals occurs. Finally, in the large distance region d ~ 5 — 6 fm, 
the a orbital becomes the fifth orbital which is the lowest one for valence nucleons around the 2a core. This situation 
is realized in the 1/2]'' state, where the two-center structure with the distance d ~ 5 fm was obtained in the present 
calculation. It is also associated with the hyperdeformation, where the declined positive-parity orbital becomes the 
fifth lowest orbital in the Nilsson's deformed shell model. In fact, the deformation /3 = 0.74 of the 1/2^ state exceeds 
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the value /3 = 0.6 for the superdeformation. On the other hand, four valence neutrons are localized in one side of the 
two-center structure. It is contrast to the molecular orbital feature of the valence proton, which is moving around the 
whole system. It seems that spatial correlation of four valence neutrons is so strong that form a P3/2 shell closure. 
As a result, totally 8 neutrons separate into two groups consisting of 2 and 6 neutrons with a weak-coupling feature. 
We here give a comment on similarity of the neutron structure between the 1/2^ and 5/2^ states. Comparing the 
neutron structure of the l/2i with that of the 5/2^, the profile of the total neutron density of the 1/2^ is similar to 
that of the 5/2i described by the neutron 2huj configuration. Therefore, we can propose an alternative interpretation 
in a mean-field picture that the l/2f state is roughly described by the 3h,uj configuration with the proton Ihu and 
the neutron 2hoj excitations. In the present calculation, we suggest that such the exotic state, the proton intruder 
state of neutron-rich nuclei, may exist as the lowest 1/2+ state at = 8 MeV. Recently, the state at 4.83 MeV 
has been assigned to be a 1/2+ state by '^He(^^Be,^^B7)X exDeriments[26|. Because of strong production via the 
proton-transfer to the "'^^Be(0+) state and analysis of angular dependence, Ota et al. suggested this 1/2+ state to be 
a proton intruder state. The present prediction of the proton intruder configuration in the 1/2+ state is consistent 
with this observation though the theoretical excitation energy of the present calculation is slightly higher than the 
experimental value. 

In the 7/2^ state in the K"^ — 1/2+ band, all the valence proton and the valence four neutrons are localized around 
one of the a cores(Fig. lUJc)). It indicates that the molecular orbital aspect disappears, while the weak-coupling 
feature of a-^Li clustering is enhanced. Because of the spatial localization, the orbitals of the valence nucleons have 
no definite parity. In this case, the unnatural parity of the total system is carried by the parity asymmetric a-^Li 
structure. It is in different situation from the 1/2+ state with the intruder proton configuration, where the parity 
originates in the proton ITkm excitation of the single-particle orbital from negative to positive one. As mentioned 
before, after the superposition of the AMD wave functions, the proton intruder state and the a+^Li cluster state 
are mixed to each other to contribute the structure change depending on J in the K'^ = 1/2+ band. The former 
component is dominant in the low spin states, while the latter component is significant in high spin states. 
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FIG. 3: Single-particle energies in the 5/2~ , 1/2+ and 7/2+ states. The energies are calculated by diagonalizing Hartree-Fock- 
hke single-particle Hamiltonian within the single AMD wave function [13] which dominates the | J^) state. 



IV. SUMMARY 



The excited states of ^■^B were studied with a method of antisymmetrized molecular dynamics (AMD). We obtained 
the largely deformed states which construct the rotational bands, K"^ — 3/2~, K"^ = 1/2+ and K"^ = 1/2". In 
these deformed states, we found various kinds of cluster aspect including molecular orbital features in the two-center 
structure. 
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FIG. 4: Distribution of the proton, neutron and matter density, and density distribution of the valence nucleons. The left three 
panels show the distribution and deformation parameter (5 of matter, proton and neutron density. In the right three panels, 
density distribution of the highest proton orbitals, the highest and the third neutron orbitals are given. The single-particle 
orbitals are extracted by diagonalizing Hartree-Fock-like single-particle Hamiltonian within a single AMD wave function [l^ . 
The results are those calculated for the intrinsic AMD wave function $amd(Z;('^), which dominates the |J^) state. The AMD 
results for the ground state of ^^Be are also shown for comparison. The box size is 10 fm. 



The K'" = 3/2^ band is the deformed one with molecular orbital structure. This band is described by the intruder 
neutron 2fiuj configuration, and interpreted as ^^Be(0]^)-|-p, where the ^^Be(0]^) has the intruder configuration and 
the additional proton strongly couples to the deformed ^^Be core. The excited neutron orbital is regarded as the a 
orbital in the molecular orbital picture. Experimentally, there is a report which has been suggested the excited state 
at 10 MeV to be a high spin state of the neutron 2hu) configuration [25], which would be a member of this K'^ = 3/2^ 
band. 

We found the proton intruder structure in the K"^ ==1/2+ band with remarkably large deformation. The deformation 
f3 = 0.74 of the band-head 1/2^ state is larger than that for the superdeformation. In the molecular orbital picture, 
the last proton is described by the a orbital. This is very exotic state with the proton intruder configuration in 
neutron-rich nuclei. In the present calculations, it was suggested that the K"^ = 1/2+ band starts from the lowest 
1/2+ state at E'x = 8 MeV. We assigned this state to the recently observed 1/2+ state at 4.83 MeV, which has been 
suggested to be the proton intruder state by Ota et al. in the '*He(^^Be,^'^B7)X experiments. In the high spin states 
of the K"^ = 1/2+ band, Q-|-^Li-like cluster structure well develops. 

It is striking that such the deformed states with highly excited configurations or well-developed cluster states appear 
in the energy region compatible to the normal excited states, even though ^^B has the neutron-shell closure in the 
ground state. Especially the molecular orbital, cr, is formed in the deformed states of ^"^B as well as neutron-rich Be 
isotopes. The cluster aspect of the deformed states in -'^'^B can be understood in natural extension of cluster structure 
of Be isotopes. That is to say, the formation of the 2a-cluster core and the role of the valence nucleons(one proton and 
four neutrons) are key in the largely deformed states of ^^B. It is challenging to investigate possible exotic structure 
with clustering in excited states of further neutron-rich B isotopes like ^^B and ^^B. 



Acknowledgments 

The authors would like to thank Dr. Ota and Prof. Shimoura for the valuable discussions. In fact, this study has 
been motivated by their suggestion of the "proton intruder state" . They are also thankful to members of Yukawa Insti- 
tute for Theoretical Physics(YITP) and Department of Physics in Kyoto University. The computational calculations 



8 



in this work were performed by the Supercomputer Projects of High Energy Accelerator Research Organization(KEK) 
and also the super computers of YITP. This work was supported by Grant-in- Aid for Scientific Research Japan So- 
ciety for the Promotion of Science and a Grant-in- Aid for Scientific Research from JSPS. It is also supported by the 
Grant-in- Aid for the 21st Century COE "Center for Diversity and Universality in Physics" from MEXT. Discussions 
during the workshop YITP-W-06-17 on Nuclear Cluster Physics held in YITP were useful to complete this work. 



[1] H. Iwasaki et al, Phys. Lett. B — bf 481, 7 (2000). 

[2] A. Navin et al, Phys. Rev. Lett. 85, 266 (2000). 

[3] S. D. Pain et al, Phys. Rev. Lett. 96, 032502 (2006). 

[4] T. Otsuka and N. Fukunishi, Phys. Rep. 264, 297 (1996). 

[5] N. Ragaki and S. Okabe, Phys. Rev. C 61, 044306 (2000); N. Ragaki, S. Okabe and K. Ikeda, Phys. Rev. C 62, 034301 

(2000). 

[6] Y. Kariada-En'yo and H. Horiuchi, Phys. Rev. C 66, 024305 (2002). 

[7] Y. Kanada-En'yo and H. Horiuchi, Phys. Rev. C 68, 014319 (2003). 

[8] W. von Oertzen, M. Freer and Y. Kanada-En'yo, Phys. Rep. 432, 43 (2006). 

[9] S. Okabe, Y. Abe and H. Tanaka, Prog. Theor. Phys. 57, 866 (1977). 

[10] M. Seya, M. Kohno, and S. Nagata, Prog. Theor. Phys. 65, 204 (1981). 

[11] W. von Oertzen, Z. Phys. A 354, 37 (1996); 357, 355(1997). 

[12] D. Scharnwebcr, U. Mosel, andW. Greiner, Phys. Rev. Lett. 24, 601 (1970). 

[13] K. Arai, Y. Ogawa, Y. Suzuki and K. Varga , Phys. Rev. C 54, 132 (1996). 

[14] A. Dote, H. Horiuchi, and Y. Kanada-En'yo, Phys. Rev. C 56, 1844 (1997) 

[15] Y. Kanada-En'yo, H. Horiuchi and A. Dote, Phys. Rev. C 60, 064304(1999). 

[16] Y. Ogawa, K. Arai, Y. Suzuki, and K. Varga, Nucl. Phys. A673 122 (2000). 

[17] M.Ito, K.Kato and K.Ikeda, Phys. Lett. B 588, 43 (2004). 

[18] M.Ito and Y.Sakuragi, Phys. Rev. C 62, 064310 (2000). 

[19] P. Descouvemont and D. Baye, Phys. Lett. B 505, 71(2001). 

[20] Y. Kanada-En'yo, H. Horiuchi and A. Ono, Phys. Rev. C 52, 628 (1995); Y. Kanada-En'yo and H. Horiuchi, Phys. Rev. 
C 52, 647 (1995). 

[21] Y. Kanada-En'yo and H. Horiuchi, Prog. Theor. Phys. SuppL142, 205 (2001). 

[22] Y. Kanada-En'yo, M. Kimura and H. Horiuchi, Comptes rendus Physique VoL4, 497 (2003). 

[23] Y. Kanada-En'yo, Phys. Rev. Lett. 81, 5291 (1998). 

[24] Y. Kanada-En'yo and H. Horiuchi Phys. Rev. C 55, 2860 (1997). 

[25] R. Kalpakchicva et al, Eur. Phys. J. A 7, 451 (2000). 

[26] Ota et al, private communication, JPS meeting, March 2007, Tokyo, Japan. 

[27] Y. Kanada-En'yo, Prog. Theor. Phys. 117, 655 (2007). 



